Experimental evidence of non-Gaussian fluctuations near a critical point. 
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The orientation fluctuations of the director of a liquid crystal are measured, by a sensitive po- 
larization interferometer, close to the Freedericksz transition, which is a second order transition 
driven by an electric field. We show that near the critical value of the field the spatially averaged 
order parameter has a generalized Gumbel distribution instead of a Gaussian one. The latter is 
recovered away from the critical point. The relevance of slow modes is pointed out. The parameter 
of generalized Gumbel is related to the effective number of degrees of freedom. 
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The fluctuations of global quantities of a system 
formed by many degrees of freedom have very often a 
Gaussian probability density function (PDF). This re- 
sult is a consequence of the central limit theorem, which 
is based on the hypothesis that the system under con- 
sideration may be decomposed into many uncorrelated 
domains. However if this hypothesis is not satisfied then 
the PDF of global quantities may take a different form. A 
few years ago it has been proposed [H, Q, 01, EL IS Q that in 
spatially extended systems, where the correlation lengths 
are of the order of the system size, the PDF, P a (x}±, °f 
a global quantity x takes under certain conditions [6( a 
form which is very well approximated by : 

a a b a 

p a(x) = -prr-T exp{-a [b a (x - s a ) - exp(-6 a (x - s a ))]}. 
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The only free parameter of P a {x) is a because b a and s a 
are fixed by the mean (x) and the variance cr^ of x '■ 
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where T(a) is the Gamma function. This distribution 
P a {x)i named the generalized Gumbel distribution (GG), 
is for a integer the PDF of the fluctuations of the a th 
largest value for an ensemble of N random and identi- 
cally distributed numbers. Instead the interpretation of a 
non integer is less clear and has been discussed in ref . pj . 
For a — 7r/2, the distribution P a is approximately the 
BHP distribution (from Bramwell, Holdsworth, Pinton). 
It has been shown in ref. [B], 0] that the GG appears in 
many different physical systems where finite size effects 
are important. An example of these non Gaussian fluc- 
tuations is the magnetization of the two dimensional XY 
model which presents a Kosterlitz-Thouless transition as 
a function of temperature. When the control parame- 
ter is close to the critical value, the correlation length 
of the system diverges and when it becomes of the order 
of the system size, then the PDF of the fluctuations of 
the magnetization has a GG form instead of the Gaus- 
sian one [H . Several other examples where the GG gives 



good fits of the PDF of the fluctuations of global param- 
eters are : the magnetization in Ising model close to the 
critical temperature, the energy dissipated in the forest 
fire model, the density of relaxing sites in granular media 
models and the power injected in a turbulent flow and 
in electroconvection P, IE IS II 13, IS 0, 11 • Except for 
the two last examples which use experimental data all of 
the other mentioned results are obtained on theoretical 
models. Therefore it is of paramount importance to check 
whether the above mentioned theoretical predictions on 
GG can be observed experimentally in other phase tran- 
sitions. We report in this letter the first experimental 
evidence that close to the critical point of a second or- 
der phase transition, the PDF of a spatially averaged 
order parameter takes the GG form when the correlation 
length is comparable to the size of the measuring region. 
The Gaussian distribution is recovered when the system 
is driven away from the critical point. We also stress that 
the deviation to the Gaussian PDF are produced by very 
slow frequencies. 

In our experiment, these properties of global variables 
have been studied using the Freedericksz transition of a 
liquid crystal (LC) submitted to an electric field E [jlficij]. 
In this system the measured global variable x is the spa- 
tially averaged alignment of the LC molecules, whose 
local direction of alignment is defined by the unit vec- 
tor ft. Let us first recall the general properties of the 
Freedericksz transition. The system under consideration 
is a LC confined between two parallel glass plates at a 
distance L (see fig. [T]). The inner surfaces of the confin- 
ing plates have transparent ITO electrodes, used to ap- 
ply the electric field Furthermore the plate surfaces, are 
coated by a thin layer of polymer mechanically rubbed in 
one direction. This surface treatment causes the align- 
ment of the LC molecules in a unique direction parallel 
to the surface (planar alignment), i.e. all the molecules 
have the same director parallel to x-axis and n — (1, 0, 0) 
(see fig. [U [Til] - The LC is submitted to an electric field 
perpendicular to the plates. To avoid the electrical po- 
larization of the LC, the electric field has a zero mean 
value which is obtained by applying a sinusoidal volt- 
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FIG. 1: a) The geometry of Freedericksz transition : director 
configuration for Vo < V c and director configuration for Vo > 
V c - b) Definition of angular displacement 9 of one nematic n. 
c) Experimental setup. A polarized laser beam is focused into 
the LC cell and a polarization interferometer measures the 
phase shift $ between the ordinary an extraordinary rays [l7| . 



age V at a frequency of 1 kHz between the ITO elec- 
trodes, i.e. V = V2Vbcos(27r ■ 1000 • t) [I EH- When 
Vo exceeds a critical value V c the planar state becomes 
unstable and the LC molecules, except those anchored 
to the glass surfaces, try to align parallel to the field, 
i.e. the director, away from the confining plates, ac- 
quires a component parallel to the applied electric field 
(z-axis) (see figOJi)). This is the Freedericksz transition 
which is a structural transformation whose properties are 
those of a second order phase transition [ij E3|- For 
Vq close to V c the motion of the director is character- 
ized by its angular displacement 9 in xz-plane (fig. Hb)), 
whose space-time dependence has the following form : 
9 = 9 (x,y,t) ain(^) @, EE El- If #0 remains small 
then its dynamics is described by a Ginzburg-Landau 
equation and one expects mean-field critical phenom- 
ena [^, E3, EH j in which 9q is the order parameter and 
v 2 

e = p% — 1 is the reduced control parameter. The global 
variable of interest is the spatially averaged alignment 

of the molecules, precisely x = \ Jo < (1 — n x) >xy 
dz ~ jj A 9 2 .dxdy/A, where A = ttD 2 /4 is the area of 
the measuring region of diameter D in the (x, y) plane 
and < . > xy stands for mean on A. As \ is a global 
variable of this system, its fluctuations, induced by the 
thermal fluctuations of 9 , depend on the ratio between 
D and the correlation length £ of 9q. The angle 9q is a 
fluctuating quantity whose correlation length and corre- 
lation time are respectively: £ = L^-^/e) -1 and r = ro/e 
where tq is a characteristic time which depends on the 
LC properties and L 2 0, EH EH • Many aspects of the 
director fluctuations, such as power spectra and correla- 
tion lengths, at Freedericksz transition have been widely 
studied both theoretically d E3, E3 and experimen- 



tally [U, El EI El- However the statistical properties 
of the spatially averaged director fluctuations have never 
been characterized as a function of the ratio N e ff = D/£. 



As this ratio is the key parameter of our study, we have 
performed the experiment in cells with three different 
thickness L = 25/jto, L = 20/ito and L = 6.7/ito. The 
results reported here are mainly those of the thinner cell 
and a detailed comparison with those of the others will 
be the aim of a longer paper. The cells are filled by a 
LC having a positive dielectric anisotropy e a (p-pentyl- 
cyanobiphenyl, 5CB, produced by Merck). For this LC 
V c = 0.720 V and r c = 55 ms in the cell with L — 6.7 ^m. 

Let us describe now how \ has been measured. The 
deformation of the director field produces an anisotropy 
of the refractive index of the LC cell. This optical 
anisotropy can be precisely estimated by measuring the 
optical path difference $ between a light beam crossing 
the cell linearly polarized along x-axis (ordinary ray) and 
another beam crossing the cell polarized along the y-axis 
(extraordinary ray). The experimental set-up employed 
is schematically shown in fig. [TJ:). The beam is pro- 
duced by a stabilized He-Ne laser (A = 632.8 nm) and 
focused into the liquid crystal cell by a converging lens 
(focal length / = 160 mm). A second lens with the same 
focal length is placed after the cell to have a confocal op- 
tical system, which insures that inside the cell the laser 
beam is parallel and has a diameter D of about 125 /im. 
The beam is normal to the cell and linearly polarized at 
45° from the x-axis, i.e., can be decomposed in an ex- 
traordinary beam and in an ordinary one. The optical 
path difference, between the ordinary and extraordinary 
beams, is measured by a very sensitive polarisation in- 
terferometer 



is given by : 



17]. After some algebra the phase shift $ 
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with (n a , n e ) the two anistotropic refractive indices [9|, 
Eo[. In term of x, we get : 
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The phase $, measured by the interferometer, is acquired 
with a resolution of 24 bits at a sampling rate of 1024 
Hz. The instrumental noise of the apparatus [I3| is three 
orders of magnitude smaller than the amplitude <5$ of 
the fluctuations of $ induced by the thermal fluctuations 
of X- 

We first check the accuracy of our experimental setup 
by measuring the time-average (x) of the global variable 
X and compare it to the results of a mean-field theory. 
In figure [2 we plot the measured (x) versus the control 
parameter e. (x) vanishes for e < and increases for 
e > 0. The experimental results are in very good agree- 
ment with theoretical predictions based on the Ginzburg- 
Landau equation using the physical properties of this LC 
without adjustable parameters. This excludes the pres- 



3 




ict 1 io u 
f (Hz) 



FIG. 2: a) Average value of the the global variable x as a 
function of e (o). Continuous line is a theoretical prediction 
based on the Ginzburg-Landau equation for 6a using the val- 
ues of this LC, with no adjustable parameters, b) Power spec- 
trum S x and the Lorentzian fit (continuous line) measured at 



e = 2 ■ 10~ 3 (o) and at e = 0.16 (□). 



ence of the weak anchoring effects described in litera- 
ture If], 11 1 ■ We observe that the model is valid even for 
large values of e. The rounding near the transition is a 
finite size effect because in our experiment N e f / — 1 for 
e ~ 0. Indeed cells with the larger L show even a more 
pronounced rounding effect close to the transition. 

To shed some light on the dynamics of the fluctuations, 
we first measure the power spectral density S x of x- As 
the slow thermal drift of the interferometer may perturb 
the statistics of the acquired signals, \ lSi high-pass fil- 
tered at 2 mHz. The power spectrum S x , measured at 
e = 0.16 and e = 0.002, are plotted in figure [H They can 
be fitted by a Lorentzian for e > : 



S x = 



So(e) 



(•5) 



^(e) represents the amplitude of fluctuations and / c (e) 
is proportional to the inverse of the relaxation time r(e) 
of #0 (fc = (7rr) _1 ). This form is the same found by 



Galatola for light-scattering measurements [14|, [15j but 
we have increased the resolution at low frequencies of 
about three orders of magnitude. The values of Sq and 
f c are obviously dependent on e and its sign. For e < 0, 
S x is the sum of two Lorentzian functions with two cut-off 
frequencies. Each frequency corresponds to a relaxation 
of the director of the LC in two different directions. The 
lowest frequency, which corresponds to 9, depends on e 
contrary to the other frequency. The cut-off frequency 




at e ~ 0.16, 8-10" 



FIG. 3: a),b),c),d) PDFofy : 
10 -3 , and 2 ■ 10~ 3 respectively. Dashed line is a Gaussian fit. 
In b), c) and d) the continuous lines are the GG distributions 
with a = 23.5, 6.6, and 2.95 respectively. 



f c decreases with e with a linear behavior as predicted 
by the Ginzburg-Landau model, i.e. 1/r = e/r , where 
the value of To agrees with that obtained from the LC 
parameters. The amplitude So has a complex dependence 
on e. This dependence, which can be understood on the 
basis of the Ginzburg-Landau model, is not relevant for 
the results presented in this letter and will be discussed 
in a longer report. 

We now turn to the main point of this letter that is 
the statistical description of the fluctuations of x- We 
consider the normalized order parameter : y = 
where a 2 is the variance of x- The probability density 
functions of y are plotted in Fig. [3] for three different 
values of e. We find that far from the critical value (e = 
0.16) the distribution is Gaussian (Fig. [5^)). In contrast, 
for a value of e closer to 0, typically e ~ 2 ■ 10~ 3 , the PDF 
of fluctuations of x are n °t Gaussian as it is clear from 
figCUd). In figure [3b) and c), we plot the distribution of 
X for two intermediate values of e. The exponential tail 
becomes more pronounced when e decreases. We want 
now to compare this distribution to a GG (eq. {!])). The 
value of the free parameter a is given by the skewness of 
the fluctuations [4[ : 

7 - (V 3 ) - -(^^)/(^^) 3/2 - -1/VE (6) 



da 4 



We obtain a = 2.95 at e - 2-10~ 3 , a = 6.6 at e - 4-10~ 3 
and a = 23.5 at e ~ 8 • 10~ 3 . Using these values in 
eq. (HJ) , we get the PDFs plotted in fig [3] as continuous 
lines which agree quite well with the experimental dis- 
tributions. The observation of the GG for e very close 
to is the main result of this letter. One may wonder 
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FIG. 4: a) The skewness 7 of the fluctuations at e = 2 ■ 10~ 3 
is plotted as a function of /hp- Inset: 7 -1 is linear in /hp//c- 
b) 7 as a function of e. Inset: 7" is linear in e. 



why the GG is observed in our experiment and not in 
other experiments on phase transitions. To answer to this 
question, let us first consider the slow modes of \ whose 
relevance for the GG distribution has been pointed out 
in ref. 5] . To confirm this point the time evolution of \ 
acquired ate = 2-10~ 3 is high-passed filtered at various 
cut-off frequencies /hp- The skewness 7 of the filtered 
signal is plotted as a function of /hp (figure Ek.)). When 
/hp is increased, we see that the skewness decreases (7 _1 
is linear in /hp). A Gaussian behavior is retrieved for 
/hp > 10/c — 0.1 Hz. These experimental results indi- 
cate that the slow modes, with frequency lower than f c , 
are responsible for the non-gaussian PDF of this global 
parameter. Previous experiments on Freedericksz transi- 
tion did not have a sufficient resolution at low frequencies 
and they erased this effect. Let us now consider the cor- 
relation length £ of #0 in the plane (x,y). This correlation 
length has to be compared with the diameter of the mea- 
suring volume, that, in our experiment, is determined by 
the laser beam diameter D inside the cell. At e = 0.002, 
we find £ = 47 /Ltm, that is £ ~ D/3. In other words 
the laser detects the fluctuations of only a few coherent 
domains and, in agreement with the theoretical predic- 
tions, these fluctuations have the GG distribution. The 
effective number of degrees of freedom of the system is 
related to the ratio N e g = D/£ oc y^e. In figure [3b), 
we plot the values of the skewness as a function of e. We 
observe that 7 goes to zero for increasing e and the Gaus- 
sian behavior is retrieved for e > 0.03. The inverse of 7 is 
linear in e, that is 7 _1 = — ^/a = p + qe = p + qN^ s . We 



measure p — —0.51 and q — —521. Thus the free param- 
eter a of the GG is a measure of the effective number of 
degrees of freedom as underlined in ref. 0, EH . For the 
magnetization of the two dimensional XY model it has 

-1 K „, R 



been found that: 7 



f 



f.f 



The dependence of \fa on N e e is the same than in our 
experiment but the coefficients depend on the system. 
As the £ is proportional to the cell thickness, we have 
verified that for cells having larger L, the GG is obtained 
for larger values of e. This is indeed the case and the 
detailed description of the results of the other cells will 
be the subject of a long article. 

In conclusion, we have experimentally shown, using the 
Freedericksz transition of a LC, that in a second order 
phase transition the fluctuations of a spatially extended 
quantity have a GG distribution if the coherence length 
is of the order of size of the measuring area. The slow 
modes, corresponding to large scales, are responsible for 
this non-Gaussian behavior. This observation confirms 
several theoretical predictions on GG, which have never 
been observed before in an experiment on a phase tran- 
sition. 
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